Study of the correlation effects in Yb + and implications for parity violation 
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Calculation of the energies, magnetic dipole hyperfine structure constants, El transition ampli- 
tudes between the low-lying states, and nuclear spin-dependent parity-nonconserving amplitudes 
for the 2 Si/2 — 2 D3/2,s/2 transitions in 171 Yb + ion is performed using two different approaches. 
First, we carried out many-body perturbation theory calculation considering Yb + as a monovalent 
system. Additional all-order calculations are carried out for selected properties. Second, we carried 
out configuration interaction calculation considering Yb as a 15-electron system and compared the 
results obtained by two methods. The accuracy of different methods is evaluated. We find that 
the monovalent description is inadequate for evaluation of some atomic properties due to significant 
mixing of the one-particle and the hole-two-particle configurations. Performing the calculation by 
such different approaches allowed us to establish the importance of various correlation effects for 
Yb + atomic properties for future improvement of theoretical precision in this complicated system. 

PACS numbers: 31. 15. am, 11.30.Er, 32.10.Fn, 32.70.Cs 
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I. INTRODUCTION 

The Yb + ions have been a subject of heightened in- 
terest in recent years owing to use of this system in a 
number of different applications including quantum in- 
formation studies 0J2], searches for variations of fun- 
damental constants |3[, and development of the optical 
frequency standards [4}tZ|- 

Manipulation and detection of a trapped Yb + hyper- 
fine qubit was described in pj. An efficient scheme to 
carry out gate operations on an array of trapped Yb + 
ions was suggested in 2]. Yb + is of particular inter- 
est to the atomic clock development due to the avail- 
ability of two different (quadrupole [j| and octupole Q) 
metastable transitions that can be used as optical fre- 
quency standards. In 2012, the performance of the op- 
tical frequency standard based on electric-octupole tran- 
sition 2 Si/ 2 (F = 0) ->• 2 Fy 2 (F = 3) in a single trapped 

Yb + ion was investigated Q- This work has demon- 
strated that the octupole transition in m Yb + can be 
used to realize an optical clock with a systematic uncer- 
tainty of 7.1 x 1CP 17 [7] ■ Moreover, it has been shown that 
a clock based on a linear combination of the quadrupole 
and the octupole transition frequencies of Yb + can have a 
significantly reduced blackbody shift 8]. An availability 
of two metastable transitions suitable for the develop- 
ment of the precision frequency standard made Yb + an 
attractive candidate for the search of the variation of the 
fine-structure constant. 

The 2 S'i/2 —> 2 Z?3/o transition in Yb + was also pro- 
posed [9j for study of the nuclear spin-dependent (SD) 
parity-nonconserving (PNC) effects. Such an experiment 
will be able to yield the nuclear anapole moment that 
arises due to parity-violating interaction between nucle- 
ons in the nucleus 10]. Study of the weak hadronic in- 



teractions is of particular interest due to significant dis- 
crepancy between constraints on weak nucleon-nucleon 
couplings obtained from the cesium anapole moment and 
those obtained from other nuclear parity-violating mea- 
surements [n], [12] . 

Accurate calculation of Yb + properties is very difficult 
owing to the large number of low-lying states of the hole- 
two-particle configurations such as 4/ 13 5d6s and their 
strong mixing with one-particle (monovalent) configura- 
tions, such as 4/ 14 6p. The properties of ytterbium ions 
were studied in a number of theory papers (see, e.g., (l3| 
and references therein). Because the main configuration 
of the ground state of Yb + is 4/ 14 6s, this ion can be con- 
sidered as a system with one electron above the closed 
shells. Alternatively, one can treat the 4/ electrons as 
the valence electrons, and consider Yb + as a system with 
15 valence electrons. Both approaches have advantages 
and drawbacks. In the following we refer to them as a 
single-electron approach and a many-electron approach. 

The advantages of a monovalent (single-electron) 
method are high accuracy and relative simplicity. In par- 
ticular, the core-valence correlations can be treated very 
accurately. However, the problem is that the states be- 
longing to the configurations with the unfilled 4/ shell, 
such as 4/ 13 6s 2 and 4/ 13 5c?6s, are lying very low in Yb + . 
A knowledge of their properties is important for a num- 
ber of experimental schemes mentioned above. A single- 
electron method is unable to treat such states since these 
are not monovalent states. Moreover, the states with 
filled 4/ shell (like 4/ 14 6p 2 Py 2 ) can strongly interact 
with a closely located state with the unfilled 4/ shell. 
This mixing can significantly affect the properties of both 
states. Again, the single-electron approach does not take 
into account this interaction that drastically affects the 
accuracy of this approach for the states where this mixing 
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is large. This effect is illustrated on the example of cal- 
culation of the magnetic dipole hyperfine structure (hfs) 
constant A of the 4/ 14 6p 2 Py 2 state. It was calculated 
by several different methods that considered Yb + as a 
monovalent system [L3l - [l6j ]. 

The resulting values are in reasonable agreement with 
each other but are factor of two smaller than the exper- 
imental result. As we show in this work, the reason for 
this discrepancy of theory and experiment is the strong 
configuration interaction between the 4/ 14 6p 2 Py 2 ancl 
4/ 13 5d6s 3 [3/2]° /2 states. 

The many-electron methods, such as the conventional 
configuration interaction (CI), allow us to study the prop- 
erties of the states with both filled and unfilled 4/ shell 
on equal footing. It also allows to take into account the 
configuration interaction between nearby levels. How- 
ever, the accuracy of the 15-electron CI approach is gen- 
erally lower than that of single-electron methods due to 
omission of the correlation corrections between the core 
[Is 2 , 5p 6 } electrons and the valence electrons. So far, 
it has not been possible to incorporate successfully the 
core-valence correlations into a many-electron CI. 

In this work, we carried out calculation of Yb + prop- 
erties using both the single-electron approach, with both 
many-body perturbation theory (MBPT) and all-order 
methods, and the 15-electron configuration interaction 
method. The use of the both approaches allows us 
to study the properties of all low-lying states. Since 
these methods are to some extent complementary to 
each other, they give us clearer picture of the impor- 
tance of various correlation effects and provide more com- 
plete theoretical description of the Yb + properties. This 
work will allow to outline a pathway for the develop- 
ment of more accurate approaches for the calculation of 
Yb + properties of interest to applications listed above. 
Because of the importance of the Yb + for various appli- 
cations, experimental values of other properties will be- 
come available in the future for further theory tests. Yb + 
is an excellent system for benchmark tests of new theo- 
retical approaches capable of describing strong electron 
correlations. A development of such new approaches is 
also needed to improve theoretical description of complex 
atoms, such as Dy or Ho, which is becoming more impor- 
tant owing to recent experimental developments and new 
proposals with these systems p7l - [l9l |. 

Another goal of this paper is to evaluate spin- 
dependent parity-violating amplitudes for the 
4/ 14 6s 2 Si /2 - 4/ 14 5d 2 D, ; transitions in Yb+ and 
to study the effects of various correlation corrections to 
this quantity. The calculation of the PNC amplitude 
is required to analyze the experimental PNC studies 
and extract the anapole moment (unless the measure- 
ments are carried out with several isotopes). Such 
experimental study with Yb + is presently underway 
in Los Alamos [9(. So far, a non-zero anapole mo- 
ment was observed only in Cs [10]. The Cs result is 
in disagreement with the nuclear physics predictions 
for the Cs anapole moment and constraints on weak 



nucleon-nucleon couplings [111 . Il2j prompting further 
investigations. The spin-dependent PNC effects in the 
4f 14 5d 2 D 3/2 - 4/ 14 6s transition of Yb+ were 

recently investigated in [l5[. The authors of [l5[ treated 
Yb + as a monovalent system. They noted significant 
cancelation between different terms contributing to the 
SD PNC amplitudes which merited further investigation 
carried out here. We note that the total uncertainty in 
the value of the anapole moment that can be extracted 
from the experiment with a single isotope includes the 
theoretical and experimental uncertainties. 

Other PNC experiments that are presently underway 
include studies with Yb [2(| , Fr [2l[ , and Ra + [22j • Large 
atomic parity violation effect was observed in neutral Yb 

The paper is organized as follows. The Sec. |H] is de- 
voted to the single-electron approach. We present the 
results of calculations of the low- lying energy levels, the 
magnetic dipole hfs constants, El transition amplitudes 
between the low-lying states, and nuclear spin-dependent 
parity-nonconserving amplitudes for the 2 Si/2 — 2 C ) 3/2,5/2 
transitions. The all-order results are also given for the en- 
ergy levels and electric dipole matrix elements. In Sec. lIIII 
we present the energy levels, magnetic dipole hfs con- 
stants, and El transition amplitudes for the low- lying 
states calculated in the framework of the 15-electron CI 
method. We compare the results obtained by either 
method. If the results differ from each other the reasons 
are analyzed. The conclusions are described in Sec. IIVI 
We use atomic units % = |e| = m e = 1 thorough the 
paper unless stated otherwise. 

II. SINGLE-ELECTRON APPROACH AND 
RESULTS 

In this approach, the 4/ electrons are considered as the 
core electrons. We start from the solution of the Dirac- 
Fock (DF) equations carrying out the self-consistency 
procedure for the [ls 2 ,...,4/ 14 ] core electrons. Then, the 
valence orbitals 6 — 8s, 6 — 7p, and 5 — 6d were constructed 
in the V^ -1 approximation (N is the total number of the 
electrons in the system). The basis set used in calcula- 
tions included virtual orbitals up to 23s, 23p, 23d, 22/, 
and 14<7 formed with the help of the recurrent procedure 
described in Refs. [H, HJ]. The MBPT corrections can 
be included by solving the equation 

H s ip n = £ n ip n (1) 

with the effective Hamiltonian defined as H c g = Hq + S, 
where Ho is the Dirac-Fock Hamiltonian and the operator 
£ takes into account virtual core excitations. 



A. Energy levels 

First, we find the energies of the low- lying states in 
various approximations and compare them with the ex- 
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TABLE I: The comparison of the energy levels calculated in 
different approximations with experiment [2f|. The ioniza- 
tion potential is given in the first line (in a.u.), the energy 
levels of the excited states are counted from the ground state 
(in cm- 1 ). The columns labeled "DF" and "MBPT" cor- 
respond to the Dirac-Fock and DF+MBPT approximations 
with the MBPT corrections included in the second order. The 
higher orders of the MBPT are included in the results listed in 
the column labeled "MBPT(HO)". The results of the single- 
double all-order calculation are listed in the column labeled 
"All-order". 





DF 


MBPT 


MBPT(HO) All-order Experiment 




0.41366 


0.45211 


0.44473 


0.45090 


0.44775" 


D3/2 


24272 


24450 


22711 


22820 


22961 


2 D 5/2 


24752 


25952 


24178 


24261 


24333 


2r>o 
-fl/2 


24702 


28636 


27945 


27945 


27062 


J 3/2 


27513 


32242 


31403 


31481 


30392 



"This is equal to the ionization potential = 98269 cm 1 [2 



perimental values [25j . As we already mentioned, we are 
able to obtain only the energies for the states with filled 
4/ shell in the framework of this approach. In Table U 
we present the ionization potential (line 1) and the en- 
ergies of the low-lying states obtained in different ap- 
proximations. At the DF stage of the calculations even 
the order of the levels is incorrect. An inclusion of the 
second-order MBPT corrections restores the correct or- 
der of levels listed in Table [I] The second-order MBPT 
values are listed in column labeled "MBPT" . 

As expected, the agreement between the experimen- 
tal and theoretical energies significantly improves with 
the inclusion of the correlation corrections beyond DF 
approximation. At the same time, the correlations are 
large and accounting for only the 2-nd order MBPT cor- 
rections is not sufficient. To calculate the energy levels 
(and subsequently other properties) more accurately, we 
need to take into account the higher orders (HO) of the 
MBPT. We designate this approximation as MBPT(HO) 
and label the results of such calculations accordingly in 
the text and the tables. 

In this approach, we include higher-order corrections 
by introducing screening coefficients Ck for the Coulomb 
lines in self-energy diagrams (see, e.g., HH). The latter 
can be calculated as an average screening of the two- 
electron Coulomb radial integrals of a given multipolar- 
ity k. These coefficients serve as an approximation to 
the insertion of polarization operator in Coulomb lines. 
The coefficients Ck were chosen as follows: Cq = 1.3, 
C\ = 0.75, and Ck — 1 for k > 2. The resulting energies 
are listed in Table U in column labeled "MBPT(HO)" . 
The ionization potential obtained in the MBPT(HO) ap- 
proximation agrees with the experiment at the level of 
0.07%; the energies of the even-parity states are within 
1% from the experimental values, and the energies of the 
states were reproduced with the 3% ac- 



curacy. 



and 2 P° /2 



As a test of the MBPT(HO) approach, we also calcu- 
lated the energy levels using the linearized single-double 
coupled-cluster method (also referred to as the all-order 
method). This method allows to include the higher-order 
correlation corrections in an ab initio way by effectively 
summing the dominant correlation contributions to all 
orders of the perturbation theory. The single-double all- 
order method was demonstrated to produce very accu- 
rate atomic properties for alkali-metal atoms and other 
monovalent systems. We refer the reader to the review 
27} for a description of the all-order approach and its 
applications. The all-order data are listed in column la- 
beled "All-order" in Table 0] These energy values have 
been previously listed in [28| ■ We find that ab initio all- 
order energy levels are close to the MBPT(HO) values 
serving as an additional verification of the MBPT(HO) 
approximation. 



B. Magnetic dipole hfs constants and El transition 
amplitudes 

To calculate magnetic dipole hfs constants and El 
transition amplitudes, we construct effective valence op- 
erators for the magnetic dipole hyperfine interaction Hm s 
and the electric dipole operator d. First, we solve the 
random-phase approximation (RPA) equations, which is 
equivalent to the summation of the dominant sequence 
of many-body diagrams to all orders of MBPT [2y, [2i| . 
Then, we include additional corrections (beyond RPA) 
to the effective operators: the core-Brueckner, structural 
radiation (SR), and normalization corrections. 

The results obtained for the hfs constants are listed 
in Table HH This table illustrates that the MBPT cor- 
rections are generally large and contribute significantly 
to the hfs constants A. The RPA, core-Brueckner, SR, 
and normalization corrections also have to be taken into 
account. The 2 Dj states are particular sensitive to dif- 
ferent corrections. For instance, the RPA correction even 
changes the sign of A( 2 D 5 /2). The SR corrections (which 
are usually relatively small) are found to be significant in 
this case and change the values of A( 2 D 3 / 2 ) an d A( 2 D 5 / 2 ) 
by more than 40%. 

The final values of all hfs constants obtained in this 
work are, in general, in reasonable agreement with the 
experimental data and other theoretical results with the 
exception of two cases. We find a significant difference 
between our value of A( 2 D 3 / 2 ) and the value found in (l5j . 
The difference is most probably due to inclusion of the 
corrections beyond RPA in the present work. Our result 
is in a good agreement with the experiment. All of the 
theoretical values are in disagreement with the experi- 
mental value of the A( 2 P°^ 2 ) demonstrating that this hfs 
constant cannot be correctly reproduced in the frame- 
work of a single-electron approach. As we will discuss 
in more detail in the section devoted to the 15-electron 
CI, this problem arises due to a strong mixing of the 
4/ 14 6p 2 P° /2 state and a nearby 4/ 14 5d6s 3 [3/2]° /2 state. 
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TABLE II: The breakdown of different contributions to the hfs constants A (in MHz) (I = 1/2, fi = 0.4919 [3fJ|). First row gives 
the DF values and the following rows give MBPT(HO), RPA, core-Brueckner (a), structural radiation (SR), and normalization 
(Norm.) corrections. The row labeled "Total" gives the final numbers. The values are compared with the experimental and 
other theoretical [Il4]l results. 







2 D 3 /2 


2 D 5/2 


2po 
/ l/2 


2po 
~3/2 


DF 


9577 


290 


111 


1542 


183 


MBPT(HO) 


2993 


109 


38 


549 


58 


RPA 


1672 


-55 


-308 


323 


132 


a 


-762 


-13 


-5 


-10 


-5 


SR 


-188 


167 


67 


-9 


-35 


Norm. 


-201 


-9 


1 


-24 


-3 


Total 


13091 


489 


-96 


2371 


330 


Experiment 


12645(2)" 


430(43) 6 


-63.6(7) c 


2104. 9(1. 3) a 


877(20) d 


Ref. [13] 


13172 






2350 


311.5 


Ref. [14] 


12730 






2317 


391 


Ref. [15] 


13217 


291 




2533 


388 


Ref. [16] 


13332 


447 


-48 


2516 


322 



"Reference O; 'Reference [3l[ 

and references therein; ^Reference [32|: ^Reference [33l| 



A possible sensitivity of A( 2 Py 2 ) to the configuration 
mixing was also mentioned in Ref. [l5j . 

We also calculated the El amplitudes for the transi- 
tions between the low-lying states and compared them 
with other available data. The lifetime t of the 2 P°n 
state was measured with a high precision in Ref. [13] to 
be equal to 8.12(2) ns. The 2 P°/ 2 state can decay to the 

2 £>3/ 2 and 2 Si/ 2 states. The decay channel to the ground 
state strongly dominates. The branching ratio from the 
2 P°/ 2 state to the metastable 2 D 3 / 2 state was measured 
to be 0.00501(15) [l[. Using two these quantities, we find 
the transition probabilities 

W{ 2 P° /2 2 S 1/2 ) = 0.995/t( 2 P° /2 ) 

= 1.23(5) x 10* s' 1 (2) 

and 

W( 2 P? /2 -> 2 D 3/2 ) = 0.00501 /t( 2 P? /2 ) 

= 6.17(18) x 10 5 s _1 . (3) 

Respectively, the "experimental" reduced matrix el- 
ements (MEs) of the electric-dipole moment oper- 
ator are K 2 Si/2||d|| 2 i\ /2 )| = 2.471(3) a.u. and 
\{*D 3/2 \\d\\*P? /2 )\ =2.97(4) a.u. 

At the present time, the most precise measurement 
of the 2 P-3/ 2 lifetime, T~( 2 Py 2 ) = 6.15(9) ns, was carried 
out in [35j ]. The 2 P^ 2 state mainly decays by the El 
transitions to the 2 Si/ 2 , 2 D 3 / 2 , and 2 D 5 / 2 states. Then 

_L_ „ Wfp* 2 Si/2) + w( 2po ^ 2 Dm) 

+ W{ 2 P° /2 2 D 5/2 ). (4) 

The transition probabilities W( 2 P°^ 2 -> 2 D 3 / 2 ) and 
W{ 2 P° /2 -> 2 D 5/2 ) were calculated in Ref. [13] to be 



3.6 x 10 5 s 1 and 1.9 x 10 5 s 1 , respectively. Thus, 
they are more than two orders of magnitude smaller than 
1/t( 2 P 3 ° /2 ) w 1.626 x 10 s s" 1 . Even if the accuracy of 
W{ 2 P° /2 -> 2 Dj) is not so high (for example, - 20-30%), 
it practically does not affect the final accuracy of the 
W( 2 Py 2 —> 2 5i/2) inferred from the experiment. Using 
the experimental value of t~( P£, 2 ), we find the probabil- 
ity of the 2 Py 2 — > 2 Si/ 2 transition from Eq. (j4]) yielding 
\{ 2 S 1/2 \\d\\ 2 P° /2 )\ fa 3.36(3) a.u. The same experimental 
value for this reduced ME was quoted in 15]. 

In Table IIII1 we present the results obtained for the 
reduced MEs of the electric dipole moment operator d in 
the DF approximation and list the MBPT(HO), RPA, 
and other corrections. We emphasize that the core- 
Bruckner, SR, and normalization corrections are small 
in this case, and we do not present them separately. The 
sum of these corrections is given in the table in the row 
labeled "Other" . We also calculated the El matrix ele- 
ments using the ab initio all-order method [27| . These 
values are listed in the row labeled "All-order" . These re- 
sults include the dominant SR, normalization, and other 
corrections to all orders. The MBPT(HO) and single- 
double all-order results are in close agreement. 

C. Parity- nonconserving amplitude 

We carried out the calculation of the spin-dependent 
PNC amplitudes for the 2 «Si/2 — > 2 ^ ) 3/2,5/2 transitions. 
The Hamiltonian describing the main part of the nuclear 
spin-dependent PNC electron-nuclear interaction can be 
written as follows 

#sd = -J=jalp[r), (5) 

where Gp w 2.2225 x 10~ 14 a.u. is the Fermi constant 
of the weak interaction, h is the dimensionless coupling 
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TABLE III: The breakdown of different contributions to the reduced MEs of the electric dipole moment operator d (in a.u.). 
First row gives the DF values. The 2nd and 3rd rows give the MBPT(HO) and RPA corrections, respectively. The row labeled 
"Other" is the sum of the core-Brueckner, structural radiation, and normalization corrections. The row labeled "Total" gives 
the final numbers. The results of the SD all-order calculation are given in the row labeled "All-order". The values are compared 
with the experimental and other theoretical [l3l [l5l . Il6l ] results. 

K 2 Sl/2]|rfi| 2 A°/ 2 )l l( 2 Sl/2lM1| 2 P 3 °/2)l l( 2 A V2 |M|| 2 A°/2)l K 2 g3/2||rf|| 2 P 3 /2)l I {^5/2 1 g l^) I 

DF 3.24 4.54 3.86 1.70 5.20 

MBPT(HO) -0.16 -0.24 -0.47 -0.23 -0.61 

RPA -0.33 -0.42 -0.32 -0.12 -0.36 

Other 0.002 -0.05 -0.006 0.001 -0.002 

Total 2.75 3.83 3.06 1.35 4.23 

All-order 2.64 3.71 2.98 1.32 

Experiment 2.471(3)" 3.36(3) 6 2.97(4)" 

Reference [13] 2.68 3.77 2.97 1.31 4.12 

Reference [15] 2.72 3.84 3.09 1.36 

Reference [16] 2.72 3.83 3.06 1.35 4.23 

"References 1. 34]; b Reference [35|. see also explanation in the text. 



constniit . o = ( ^ q J and 75 arc the Dirac matrices, 

J" is the nuclear spin, and p(r) is the nuclear density dis- 
tribution. 



We consider the nucleus to be a uniformly charged 
sphere. Then, 



p(r) = 



d(R-r). 



The root-mean-square charge radius is r rms = 5.2891 
fm [IHl and, respectively, the nuclear radius R = 
^573 r rms « 6.828 fm. 

If \i) and |/) are the initial and final atomic states of 
the same nominal parity, then taking into account the nu- 
clear SD part of the PNC interaction in the lowest nonva- 
nishing order, one can write the electric dipole transition 
matrix element as 



(/K,Sd|«) = 



E 



(f\d q \n)(n\H SD \i) 



Ei — E„ 
(f\H SJ3 \n)(n\d q \i) 



(0) 



where \a) = \J a F a M a ), F = I + J is the total angular 
momentum, M is the projection of F, and i?sD is given 
by Eq. ©. 

The expression for the reduced ME of dsn was derived 



in [37] and is given by 
(Jf F f HdsDHW 



yjl(l + 1)(27+ 1)(2F 4 + l)(2f> + 1) 

E \^y Jf - J ' { J i J i I 

(J f \\d\\n, J n )(n, J n \\H SD \\Ji) 



Jn Jf 1 

F f Fi I 



En. — Ei 



+ (-1) 



Fr-Fi 





Jf 








{t 


I 


Ff } 


\ Fi Ff 





(J f\\H SI) \\n, J n )(n, J n 



E, n - E f 



(J) 



For the 2 Si/ 2 — > 2 Dj transitions, where J= 3/2 and 5/2, 
in 171 Yb (I = 1/2) we obtain from Eq. Q 

( 2 Dj, F f \\d SD \\ 2 S 1/2 ,F t ) 



3(2F l + l)(2F / + l) 



Jn J 1 

Ff F I 



< 2 £>j||d||n,Jn)(n,Jn||gsD|| 2 gl/2) 
E n ~ E2 Sl/2 

( U F f -Fi j Jn J 1 \ / Jn 1/2 1 

+ ( j I 1/2 1/2 Fj / \ Fi F f I 
( 2 Dj\\H SD \\n,J n )(n,J n \\d\\ 2 S 1/2 ) 



(8) 



E n — E2 Dj 

For subsequent calculations it is convenient to write 

( 2 Dj, Ff ||d S D|| 2 Sl/2,^) 

= ( 2 Dj, F f \\d ■ Ri ■ H SD \\ 2 S 1/2 , F^ 

+ ( 2 D J ,F f \\H su -R 2 -d\\ 2 S 1/2 ,F l ), (9) 



6 



TABLE IV: The nuclear spin-dependent PNC amplitude (in units ix ■ 10 12 |e|oo), where ao is the Bohr radius. The values 
obtained in the DF and DF+MBPT(HO) approximations are listed in the columns labeled "DF" and "MBPT(HO)" . The 
RPA and other corrections are included in the results listed in the column labeled "RPA+other". The rows labeled "core" 
show contributions of the core excitations. The final values (given in the rows labeled "Total") are compared with the results 
obtained in Ref. [l5l ]. 



Ff 


Fi 






DF 


MBPT(HO) 


RPA+other 


Ref. [15] 


1 





( 2 D 3/ 2,Ff\\d ■ i?i ■ H SD \ 

\ Of At J | | A iJU | 


\ 2 S 1/2 ,Fi) 


6.2 


6.6 


6.9 








( 2 fl,„ Ft lli/sr) • R-y ■ d\ 


\ 2 Smo F) 








-5.1 










core 


0.7 


0.7 


0.8 








Total: ( 3 D 3/a ,F f \\dsT>\ 


I 2 5i/ 2 ,K> 


6.9 


7.3 


2.6 


3.1(1.9) 


1 


1 


{ 2 D 3 / 2 ,Ff\\d • i?i • H su \ 


| 2 5i/ 2 ,^> 


1.5 


1.6 


1.5 








( 2 D-no Ff\\H<^r, ■ Ro ■ d\ 


l 2 Si io Fi) 

^1/2 ) 1/ 








-3.2 










core 


0.2 


0.2 


0.2 








Total: ( 2 D 3/2 ,F f \\d SD \ 


| 2 Sl/2,^> 


1.7 


1.8 


-1.5 


-1.3(4) 


2 


1 


( 2 C 3 /2,f/||d-i?i -Hsnl 


| 2 5i/ 2 ,^> 


-3.3 


-3.5 


-3.6 








( 2 D 3/2 ,F f \\Hs-D-R2-d\ 


I 2 5i /2 ,K> 








1.8 










core 


-0.4 


-0.4 


-0.4 








Total: ( 2 D 3/a ,F f \\dsT>\ 


| 2 5i/ 2 ,^> 


-3.7 


-3.9 


-2.2 


-2.6(1.3) 


2 


1 


{ 2 D B/2 ,F f \\d-Ri-H so \ 


| 2 5i/ 2 ,^> 






0.2 








( 2 D 5/2 ,F f \\H B u-R 2 -d\ 


I 2 5i /2 ,K> 






-1.1 










core 






-0.1 








Total: ( 2 D 5/2 ,Ff\\d SD \ 


| 2 Sl/2,-Fi) 






-1.0 





where we denote the terms involving summations over n 
by Ri and i?2- 

To calculate the nuclear spin-dependent PNC ampli- 
tude defined by Eq. ([5]), one needs to sum over all pos- 
sible intermediate states or to solve the corresponding 
inhomogeneous equation. Here, we solve the inhomo- 
geneous equation using the Sternheimer-Dalgarno-Lewis 
method |3g, |39| in the valence sector. The results ob- 
tained in different approximations are presented in Ta- 
bleHVl 

We carried out the calculations in the DF and 
DF+MBPT(HO) (i.e., including the higher orders of the 
MBPT) approximations. Note that in these approxima- 
tions ^DjIliJsDll 2 ^,/') = and, respectively, the second 
term in Eq. (J9j) is also zero. 

Then, we solved the RPA equations which is equiva- 
lent to the summation of the corresponding many-body 
diagrams to all orders for both d and H^b operators 
in Eq. ©. Smaller contributions that include core- 
Brueckner, structural radiation, and normalization cor- 
rections were also taken into account. When the RPA 
corrections are included, the intermediate nP^ 2 states 
also contribute to the spin-dependent PNC amplitude 
drastically increasing (in absolute value) the second term 
in Eq. ©. 

The initial and final states are the many-electron 
states. Therefore, we need to account for the core excita- 
tions. This contribution (labeled as "core" in Table llV|) 
was calculated in the DF and RPA approximations. 

Table M\ (see the column labeled "RPA+other") il- 
lustrates that the two terms in Eq. © are comparable 



in their magnitude but have the opposite sign for all 
Fi -+ Ff transitions. Therefore, they partially cancel 
each other. Unfortunately, the accuracy of the calcula- 
tion of the second term is expected to be rather poor 
since the intermediate 2 Py 2 state contributes to the sec- 
ond term at the level of 90%. The quality of the wave 
function for this state near the nucleus is expected to be 
low because of the discrepancies of the theoretical and 
experimental values for the magnetic dipole hfs constant 
A( 2 P^ 2 ) discussed above. Because of the significant can- 
cellation between terms, the final numbers in Table IIVI 
are expected to only give an order of magnitude estimates 
of the spin-dependent PNC amplitudes. 

A similar single-electron approach was used by Dzuba 
and Flambaum in 15] for calculating the PNC ampli- 
tude. They have rescaled the ab initio value of the ME 
(n\H SB \m) as 

(n|i?SD|m)rescaled = \ ^ X P^| ^cxp(™) ( n [ff gD | m ) 

y Ath{n)A th (m) 

(10) 

where A exp (k) and A t h{k) are the experimental and the- 
oretical values of the magnetic dipole hfs constant of the 
state k. The assumption that 

( 2 D 3/2 \H SB \ 2 P° /2 ) 

rescaled 

~ J A exp ( 2 D 3 / 2 ) A OX p( 2 P 3 /2 ) 

may not hold for the Yb + ions due to the mixing of the 
4/ 14 6p 2 P° /2 and 4/ 13 5d6s 3 [3/2]° /2 states. An admix- 
ture of the configuration 4/ 13 5c?6s to the leading config- 
uration 4/ 14 6p of the 2 Py 2 state leads to an additional 
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contribution to the hfs constant A( 2 Py 2 ) which is pro- 
portional to (4/ 13 5d6s|i?hfs|4/ 13 5d6s). This is a large 
contribution. However, the configuration 4/ 13 5cZ6s does 
not contribute explicitly to (4/ 14 5d|iJsD|4/ 13 5d6s) be- 
cause the one-electron ME (4/|_ffsD|6s) = 0. Our values 
are in agreement with the results of Dzuba and Flam- 
baum 115] within the estimated uncertainties. 



III. 15-ELECTRON CONFIGURATION 
INTERACTION 

We demonstrated in the preceding sections that the 
single-electron method sometimes fails to correctly pre- 
dict certain properties of the Yb + ions due to mixing 
of configurations outside of the monovalent states space. 
This mixing can be taken into account within the frame- 
work of the 15-electron CI. In this approach, the 4/ elec- 
trons are also considered as the valence electrons. 

We again start from the solution of the Dirac-Fock 
equations, but the construction of the DF orbitals is 
more complicated than in the monovalent approxima- 
tion described in the preceding section. The odd-parity 
low-lying levels belong to three different configurations, 
4/ 13 6s 2 , 4/ 14 6p, and 4/ 13 5d6s. Therefore, if we con- 
struct the basis set in a standard way, i.e., in the V w_1 
approximation, the 4/ 13 6s 2 and 4/ 13 5e?6s states will have 
much higher energy than the 4/ 1 6p states and, respec- 
tively, there will be no mixing interaction between these 
configurations. To avoid this problem, we carry out the 
initial self-consistency procedure for the [ls 2 ,...,4/ 14 , 6p] 
configuration. Then, all electrons were frozen and two 
electrons (one from the 4/ shell and another one from the 
6p shell) were moved to the 6s shell. Thus, the 6s orbital 
was constructed for the 4/ 13 6s 2 configuration. Next, all 
electrons were frozen again and one electron from the 6s 
shell was moved to the 5d shell. The 5d 3 / 2 ,5/2 orbitals 
were constructed for the 4/ 13 5<i6s configuration. 

The basis set used in the CI calculations included vir- 
tual orbitals up to 8s, 8p, 7d, 7/, and 5g. The virtual 
orbitals were constructed as described in [2J, [24| ■ As a re- 
sult, the basis set we used for these calculations is rather 
short since the size of the configuration space grows very 
rapidly with the increase of the basis set. The configu- 
ration space was formed by allowing single and double 
excitations for the even-parity states from the configu- 
rations 4/ 14 6s and 4/ 14 5<i and for the odd-parity states 
from the configurations 4/ 14 6p, 4/ 13 6s 2 , and 4/ 13 5e?6s. 
To check convergence of the CI, we calculated the low- 
lying energy levels for three cases: 1) including the single 
and double excitations to the shells 6s, Qp, 5d, and 5/ 
(we designate it [6sp5df]), 2) including the single and 
double excitations to [7sp6df5g], and 3) including the 
single and double excitations to [8sp7df5g]. In the last 
case the configuration space consisted of ~ 2 300 000 de- 
terminants and calculation of the energy levels was rather 
lengthy. 



A. Energy levels 

The low-lying energy levels were calculated using the 
three CI spaces described above. The results are pre- 
sented in Table El 

We were able to reproduce the low-lying even- and 
odd-parity states belonging to five different configura- 
tions 4/ 14 6s, 4/ 14 5d, 4/ 14 6p, 4/ 13 6s 2 , and 4/ 13 5d6s (the 
column \8sp7dfhg] in Table FVT) reasonably well. The theo- 
retical energy levels for the 4/ 14 6p 2 PJ states are located 
deeper than the experimental levels. It is not surpris- 
ing since the initial self-consistency Dirac-Fock procedure 
was carried out for this configuration. The levels of the 
Af 13 5d6s configuration are in reasonable agreement with 
the experimental data. The 4/ 14 5d 2 Dj states are lying 
5-7% higher than the experimental levels. This is also 
expected because, as we mentioned above, the 5c? orbital 
was constructed not for the 4/ 14 5d configuration but for 
the 4/ 13 5<i6s configuration. We observe the worst agree- 
ment with the experiment for the 4/ 13 6s 2 2 F!^ 2 state. A 
reason is a particular sensitivity of this state to the con- 
figuration interaction. It was confirmed by calculations 
carried out with other (smaller) sets of configurations 
(not included in Table IViT]) . We assume that more con- 
figurations have to be taken into account to reproduce 
this energy level with good precision. 

It is worth noting that it was essential to include the 
5g shell into consideration as illustrated by the compar- 
ison of the columns [6sp5df] and [7sp6df5g] in Table [V] 
Most of the observed differences in the energy levels listed 
in these two columns are due to including configurations 
involving the 5g orbitals into the CI space. A number 
of levels are very sensitive to these configurations. An 
addition of the 8s, 8p, 7d, and 7/ shells (compare the 
columns [7spQdf5g] and [8sp7dfbg] in the table) led to 
much smaller changes in the energy levels. Comparison of 
these three sets appears to indicate that further extension 
of the CI space (which will be extremely time-consuming) 
will not lead to any qualitative changes for a majority 
of the states. The results obtained in the framework of 
the single-double all-order approach are presented in Ta- 
ble IVIII for comparison. 

In Sec. im we discussed a poor agreement between the 
experimental value of the magnetic dipole hfs constant 
A(4f 14 6p 2 P£/ 2 ) an d the value obtained in the single- 
electron approach. A strong interaction of this state with 
the closely lying 4/ 13 5d6s 3 [3/2]gy 2 state was suggested 
as a possible reason of this disagreement. Our calcu- 
lation with the [8sp7df 5g] CI space reproduced the dif- 
ference between the energies of the and 3 [3/2]gy 2 
states almost perfectly (1649 cm -1 ), though the order of 
the levels was not correct. The experimental difference 
is 1634 cm . It makes us confident that the configu- 
ration mixing of these two states is taken into account 
sufficiently correctly. We would like to stress that our 
calculations are purely ab initio. No semi-empirical pa- 
rameters were used in the framework of this approach. 
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TABLE V: The energy levels of the low- lying excited states counted from the ground state (in cm -1 ). The columns [6sp5df], 
[7sp6df5g], and [8sp7df5g] give results obtained using different sets of the configurations described in the text. The results 
of the single-double all-order calculations are presented in the column labeled "All-order" for comparison. The experimental 
energy levels [IfJ are presented in the column labeled "Experiment" . 



Config. 


Term 


[6s P 5df] 


[7s P 6df5g] 


[8s P 7df5g] 


All-order 


Experiment 


4/ 14 6s 


2 Sl/2 

















4/ 14 5d 




29978 


24237 


24615 


22820 


22961 




2 D 5/2 


30283 


25068 


25464 


24261 


24333 



4/ 13 6s 2 
4/ 13 5d6s 
4/ 14 6p 
4/ 13 5d6s 

4/ 14 6p 
4/ 13 5d6s 



-^7/2 

3 [3/2]^ /2 

2po 
*l/2 

3 [9/2]^ /2 

2po 
*3/2 

3 [H/2] 11/2 
3 [H/2] 13 /2 
3 [5/2]? /2 
3 [5/2]^ /2 



24621 
22977 
21266 
26232 
27595 
24288 
27732 
28297 
29928 
30323 



26735 
25992 
24057 
28782 
30169 
27093 
30412 
31165 
32329 
32730 



26760 
26201 
24289 
28973 
30364 
27324 
30616 
31407 
32531 
32939 



27945 



31481 



21419 
26759 
27062 
28758 
30224 
30392 
30563 
31632 
31979 
32731 



TABLE VI: The magnetic dipole hfs constants A (in MHz). The columns [6sp5df], [7sp6df5g], and [8sp7df5g] give results 
obtained using different sets of the configurations described in the text. The MBPT(HO) results are presented in the column 
labeled "MBPT(HO)" for comparison. The available experimental values are given in column labeled "Experiment". 

Configuration Term [Qsp5df] [7sp6df5g] [8sp7df5g] MBPT(HO) Experiment 



4/ 14 6s 
4/ 14 5d 

4/ 13 6s 2 
4/ 13 5rf6s 
4/ 14 6p 
4/ 13 5d6s 

4/ 14 6 P 
4/ 13 5d6s 



/2 
3/2 
5/2 



r 7/2 

3 [3/2] 



2pa 



5/2 



3 [3/2]? /2 
3 [9/2]^ /2 

2po 
^3/2 

3 [H/2]? 1/2 
3 [H/2] 13 /2 
3 [5/2]? /2 
3 [5/2]? /2 



18430 
690 
252 

946 
4520 
1264 
-964 
-719 

783 
1427 
2036 
3208 
1518 



12916 
425 
154 

973 
3848 
1437 
-742 
-436 

763 
1347 
1782 
2776 
1246 



12679 
455 
164 

977 
3841 
1532 
-798 
-430 

765 
1365 
1776 
2770 
1237 



13091 

489 
-96 



2371 



330 



12645(2)° 
430(43) 6 
-63.6(7) c 

905.0(5) d 
2104.9(1.3)" 

877(20) e 



"Reference [141 ]; Reference [31[ and references therein; c Reference [32|; Reference [40|; e Reference [331 ]. 



TABLE VII: The absolute values of the reduced MEs 
I (7' I Ml 17) I ( m a - u -)> where 7 are the even-parity states and 
7' are the odd-parity states. 



/ 14 6s 2 S 1/2 f 14 5 d 2 D. 



3/2 



/ 14 5d 2 D 



5/2 



/2 
2po 
^3/2 

/ 13 6s 2 2 F 7 ° /2 
/ 13 5d6s 3 [3/2]° 
/ 13 5d6s 3 [3/2]° 



/2 
3/2 



2.51 
3.32 



1.10 



2.53 
1.05 

0.00075 
0.27 



3.27 
0.063 
0.0037 

0.86 



In the next subsection, we present the values of the 
magnetic dipole hfs constants and El transition ampli- 



tudes between the low- lying states. We compare these 
results with those obtained in the single-electron approx- 
imation and discuss the role of the mixing of monovalent 
and one-hole-two-particle configurations. 



hfs constants, El transition amplitudes, and 
other observables 



The values of the magnetic dipole hfs constants ob- 
tained using three sets of configurations are listed in Ta- 
ble IVI1 We compare these results with the values listed 
m Table [TT] obtained using the single-electron method. 
We will discuss the results obtained for the largest 
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TABLE VIII: Comparison of the reduced MEs of the electric dipole moment found in the single-electron approach (the row 
DF+MBPT) and in the 15-electron CI approach (the row 15-el. CI). The values in the row DF+MBPT include the MBPT(HO), 
RPA and smaller corrections. The experimental values are presented in third row. 

K 2 S 1/2 ||rf|| 2 P° /2 )| |( 2 Si/ 2 |M|| 2 P 3 ° /2 )| \( 2 D 3 /2\\d\\ 2 P? /2 )\ \( 2 D 3/2 \\d\\ 2 P° /2 )\ |( 2 J 5/2 ||rf|| 2 P° /2 )| 
DF+MBPT 2.75 3.83 3.06 1.35 4.23 
15-el. CI 2.51 3.32 2.53 1.05 3.27 
Experiment 2.471(3)" 3.36(3)' 2.97(4)° 

"References 1, 34]; ^Reference 35], see also explanation in the text. 



[8sp7df5g] CI space. 

For the even-parity states, the hfs constants found in 
the 15-electron CI are close to the values obtained at 
the MBPT(HO) stage (see Table HO . Such an agreement 
looks reasonable. The CI results include the correlation 
corrections between the 4/ and other valence electrons. 
In the single-electron approach these are core- valence cor- 
relations. At the same time the core excitations from 
all shells up to 4/ are completely disregarded in the 15- 
electron CI approach. The configuration mixing does 
not play very significant role for the even-parity states 
considered here because the even states with unfilled 4/ 
shell are located rather far from the ground- and 2 D 3 / 2 ,5/2 
states. 

The value of A( 2 Py 2 ) is very close to the result ob- 
tained in the single-electron approach in the DF approx- 
imation. Our analysis shows that in the many-electron 
case the contributions of all electrons (except the 6pi/ 2 ) 
nearly cancels with each other and the final value is de- 
termined almost completely by the contribution of the 
6pi/2 electron. 

The most significant disagreement between the ex- 
perimental and theoretical hfs constants in the single- 
electron approach (by a factor of 2.7) was found for 
the 2 P%i2 state. This problem is resolved in the many- 
electron calculation. The 15-electron CI gives the value 
765 MHz which differs from the experiment by only 13% 
(see Table IVll column [8sp7df5g]). The admixture of the 
4/ 13 5d6s 3 [3/2]° /2 state to the 4/ 14 6p 2 P° /2 state leads to 
an appearance of the contribution from the one-electron 
ME (6s|.Hhfs|6s) which strongly affects the value of this 
constant. Based on the results obtained in the single- 
electron approach (see Table HU the row RPA) , we es- 
timate that the core- valence correlation corrections will 
increase this number making it even closer to the experi- 
mental result. Thus, if the interaction between the 2 P£/ 2 
and 3 [3/2]°^ 2 states is taken into account, the value of 
A( 2 Py^) turns out to be in a good agreement with the 
experimental result. 

We found only one experimental result for the states 
with the unfilled 4/ shell listed in Table ED (for the 
4/ 13 6s 2 2 Fy 2 state). Our value agrees with the exper- 
iment at the level of 7%. An assignment of uncertainties 
to the hfs constants of these states is not trivial. One 
source of errors is the insufficiently large CI space. The 
corresponding uncertainties may be estimated by a com- 



parison of the results obtained for the \7 sp&dfbg] and 
[8sp7df5g] CI spaces. We see that the difference is not so 
large (at the level of a few percent). Another source of un- 
certainties is the core- valence correlations omitted in this 
approach. A magnitude of these corrections can be esti- 
mated using the results obtained in DF+MBPT method 
(see Table ILL]) . For the large hfs constants A( 2 Si/ 2 ) and 
A( 2 P° /2 ), they contribute - 30-35%. In the 15-electron 
approach these corrections are expected to be smaller be- 
cause the CI core [ls 2 ,...,5p 6 ] contains less electrons and 
is more "hard" than the [ls 2 ,...,5p 6 , 4/ 14 ] core. Thus, 
we estimate accuracy of our results to be at the level of 
25-30%. For comparison we also present in Table [VTl the 
results obtained in the framework of the single-electron 
MBPT(HO) approach. 

We also calculated the El transition amplitudes for the 
low-lying states. The reduced MEs of the electric dipole 
operator d are given in Table I VIII It is instructive to 
compare the results obtained by the single-electron and 
15-electron CI methods. This comparison is carried out 
in Table EnH 

The results indicate the following trend. The values of 
the MEs of the transitions between the ground and 2 Pj 
states are closer to the experimental values in the 15- 
electron CI approach, while DF+MBPT method gives 
better agreement with the experiment for the P°i 2 ~ 
2 Z?3/2 transition. As we noted above, taking into ac- 
count the configuration mixing is important for the 2 P$/ 2 
state. This mixing also manifests itself for the reduced 
ME |( 2 S , i/2lMI| 2 ^ , 3 / 2 )l though its influence is weaker than 
for the hfs constant A( 2 P^ 2 ). 

We conclude that all \( 2 Dj\\d\\ 2 P],}\ matrix elements 
are obtained to better accuracy in the single-electron 
method. In the single-electron approach, the 5d orbital 
was constructed for the f 14 5d configuration which is "na- 
tive" for the 2 Dj states. In the 15-electron CI approach, 
it was constructed for the f 13 5d6s configuration. Most 
likely, the set of configurations used to form the wave 
function of the 2 Dj states is not sufficiently large (even 
for the biggest CI space that we have considered) to cor- 
rectly reproduce their properties. 

In the recent work of Huntemann et al. 7], the 
quadrupole moment O of the 4/ 13 6s 2 2 Fy 2 state was 
measured to be O = —0.041(5) ea\. The quadrupole mo- 
ment O of a state \yJ) (where 7 designates all quantum 
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TABLE IX: The quadrupole moment 6 of the 4/ 13 6s 2 2 F 7 ° /2 
state (in ea ). The results are presented for three CI spaces 
[6sp5df], [7sp6df5g], and [8sp7df5g] and compared with the 
experimental and another theoretical value. 
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[6sp5df] -0.40 
[7sp6df5g] -0.06 
[8sp7df5g] -0.20 
Other theory [41] -0.22 
Experiment [7J -0.041(5) 



because the matrix element of the SD PNC Hamiltonian 
( 2 D j\\HsD\\ 2 Py-)) is greatly increased when we include 
the RPA and other corrections. To perform similar cal- 
culation in the framework of the 15-electron CI, we need 
to include the subtraction diagrams into consideration, 
what makes this approach very unstable. 

Formulating Cl+all-order approach that can treat two- 
particle-one-hole states on the same footing as the mono- 
valent states appears to be a promising way for a devel- 
opment of methodologies capable to further improve the 
calculation accuracy of the Yb + properties. 



numbers except J) is determined as 

e ^V / w^5^i^ hJ||<33lllJ> ' (12) 

where Qi is the electric quadrupole operator. 

We carried out the calculation of 6(4/ 13 6s 2 2 Fy 2 ) 
for three increasing CI spaces [6sp5df], [7sp6df5g], and 
[8sp7df5g] . The results are presented in Table IIXI As il- 
lustrated by the table, our result obtained for the largest 
CI space coincide with the theoretical value of Ref. [4J| 
and is 5 times greater (in absolute value) than the ex- 
perimental result. At the same time, we see that Q 
is very sensitive to the configuration interaction. The 
quadrupole moment is rather small due to large cance- 
lations of one-electron contributions which is expected 
to make its accurate calculation more difficult. All this 
makes the result obtained even for the biggest CI space 
[8sp7df5g] rather inconclusive. Based on our calcula- 
tions, we can only roughly estimate this quantity as 
9(4/ 13 6s 2 2 F° /2 ) ~ -0.1 eag. 

Finally, we note that an attempt to take into account 
core- valence correlations by combining the 15-electron CI 
with the MBPT was unsuccessful. The main problem 
which was repeatedly discussed earlier (see, e.g., [421 ]) 
is instability of the MBPT for the mean-field potential 
V N , which includes a large number of valence electrons. 
An accounting for the MBPT corrections leads to an ap- 
pearance of huge contribution from the subtraction di- 
agrams (43|. These diagrams are calculated only in the 
second order of the MBPT. This is insufficient for accu- 
rate treatment of the core- valence correlations. 

This problem does not allow us to calculate the SD 
PNC amplitude more accurately by the 15-electron CI 
method than it was done in the single-electron approach 



IV. CONCLUSION 

To conclude, we calculated the energies, magnetic 
dipole hfs constants, El transition amplitudes between 
the low-lying states, and the nuclear spin-dependent 
parity-nonconserving amplitudes for the 2 Si/2 — ^3/2,5/2 
transitions. Our calculations were carried out in the 
framework of the single-electron DF+MBPT method and 
by the 15-electron CI method. All-order calculations 
were also carried out for selected properties using the 
linearized single-double coupled-cluster method. 

The specific character of Yb + ion manifests itself due 
to the presence of the low-lying states with unfilled 4/ 
shell. A configuration interaction between them and the 
states with filled 4/ shell significantly affects the prop- 
erties of both types of states. We demonstrated this 
configuration mixing by analyzing the properties of the 
4/ 14 6p 2 P$/2 state. In particular, we found that an ad- 
mixture of the nearby 4/ 13 5ei6s 3 [3/2]g^ 2 state should be 
taken into account. 

Various contributions to the spin-dependent parity- 
violating amplitude are discussed and a method to im- 
prove accuracy further is proposed. 
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